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AXISYMMETRIC VIBRATIONS OF PARTIALLY 
LIQUI D-FILLED CYLINDRI CAL CONTA I NERS 

SECTION I. INTRODUCTION 


The longitudinal vibration of launch vehicles represents a significant 
problem. During ignition, lift-off or cutoff, longitudinal vibrations occur as 
transients and may cause considerable dynamic loads. However, longitudinal 
vibrations can also be excited during power flight as has been realized in the 
case of the Titan missile. Coupling of vehicle structure oscillations with pres- 
sure oscillations of the propellant system generated thrust vibrations which 
intensified the original structure vibrations. This created unstable oscillations 
of the entire system. So, attention must be paid to the dynamic behavior of 
launch vehicles in the axial direction. In particular the dynamic behavior of the 
liquid-filled flexible shell container represents a significant problem. The liq- 
uid propellant inside the tanks constitutes a high percentage of the vehicle 
masses. These masses supported by the elastic container walls act like spring 
mass systems and- -coupled- -may generate the fundamental frequencies of the 
entire system. Hence, without a clear understanding of the liquid shell inter- 
action, no exact vibration analysis of the vehicle can be performed. On the 
other hand, the pressure waves generated by this interaction may be adverse 
to propellant flow, pump performance and finally, as already mentioned, to the 
engine performance. Thus, the vibration analysis of partially liquid-filled 
shell containers is presently considered a problem of particular interest. This 
interest stimulated considerable scientific work which is cited in reference form 
[1] and specifically in other publications [2-7] . v 

Although the behavior of liquid propellant in flexible containers has been 
regarded important and considerable theoretical work already has been accom- 
plished, various problems are still left in this area and should be investigated. 
To discuss, at least, some of these problems, analytical methods have been 
used in this report based upon a very simple physical model, which is a flexible 
cylindrical container with a flat rigid bottom. Although this model appears to 
be rather unrealistic and academic it provides a good understanding of the vibra- 
tion mechanism and is convenient for computational purposes. The analyses of 
this paper are restricted to axisymmetric vibrations because of their impor- 
tance for space vehicle systems. Furthermore, the analyses are based on 
membrane theory, since- -as it is well known — axial bending of the cylinder wall 
affects only slightly the fundamental frequencies of the system. Comparison of 

Kana, D. D. ; Glaser, R. F. ; Eulitz , W. R. ; and Abramson, H. N. : Longitu- 
dinal Vibration of Spring Supported Cylindrical Membrane Shells Containing 

Liquid, (in preparation) 



experimental and theoretical data as shown in Section II confirms this behavior. 
See also Kana, Glaser, Eulitz and Abramson (see footnote on page i) and Kana 

[5]. 


Now the special problems to be considered in this report will be dis- 
cussed. First of all, emphasis is placed on the question of whether or not the 
effect of the liquid surface oscillations is significant. In other words, is it 
necessary to consider the combined liquid- shell -surface wave system or only 
the simpler liquid-shell system assuming zero pressure at the free surface? 

The modal analysis — applied for both cases (Sections in through VI) — shows 
that the amount of coupling among the "shell deflection modes" and the "liquid 
surface modes" is small. This behavior is demonstrated by the numerical 
example of the Saturn V vehicle (Section II) . In Hwang [7] the same behavior 
is shown for the liquid -filled hemispherical shell. Consequently the analysis 
based on zero ullage pressure can be considered adequate for predicting natu- 
ral frequencies and modes. This is also confirmed by the comparison of theo- 
retical and experimental results in Section n. In view of the theory of eigen- 
value problems , however , the coupled analysis represents an interesting example 
of a problem having solutions which consist of pairs of eigenfunctions. Accord- 
ingly the orthogonality relations (Section IV) appear somewhat unusual. 

Sections VI and VII discuss the influence of liquid compressibility. Al- 
though it is recognized that for large flexible containers the compressibility 
can be neglected it is of some interest to know the conditions under which com- 
pressibility becomes effective. This problem is closely related to the theory of 
Waterhammer. The third table in Section II considers experimental and theo- 
retical values of the first frequency of a liquid-filled flexible steel tube at dif- 
ferent filling heights. The theoretical values are calculated for compressible 
and incompressible liquids. In this way the effect of liquid compressibility is 
realized. Section VII contains a forced vibration analysis assuming longitudinal 
excitation of the container. It is shown that the results of this analysis are in 
agreement with that of Bleich [ 2] . 

The above mentioned experiments and many others have been performed 
by SWRI, San Antonio, Texas. The experiments have been carried out using 
cylindrical tanks of different materials, unstiffened and stiffened by rings, 
having flat rigid bottoms and also flexible flat and elliptical bottoms. Emphasis 
was placed on the axisymmetric vibrations because of their importance to launch 
vehicle vibrations. The experimental outcomes agree with those of Palmer and 
Asher [6] , i.e. , axisymmetric modes are very difficult to determine as they 
are largely obscured by the strong presence of nonaxisymmetrie response. 
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Furthermore the experimental data show that the nonaxisymmetric response 
originates from dynamic instability effects of the tank wall. A linear axisym- 
metrie response has been obtained only if considerably stiff containers as tube 
like shells, brass shells, and ring stiffened shells of other materials have been 
utilized. In this note, the experimental data for the unstiffened steel and brass 
shells having flat rigid bottoms are compared with the analytical result using 
the membrane model discussed above. 


SECTION II. NUMERICAL RESULTS, COMPARISON 
W ITH TEST RESULTS 


The analytical methods applied in this section are programmed on IBM 
709-4. However, some of the numerical results are obtained using approxima- 
tion formula as shown in Sections V and VI. 


Large Aluminum Tank 


The analysis applied is that of Section IV based on the assumption of 
coupled liquid -shell -surface oscillations. It results in two sets of eigensolu- 
tions with separated frequency ranges. The k^ 1 normal mode consists _of a 
pair of functions: shell deformation w^(z) and liquid surface elevation ^(r). 
The frequencies of the lower range correspond to modes with predominant sur- 
face elevation while the higher frequencies belong to "shell deformation modes. " 
This behavior is demonstrated by Figures 1 through 6 which show the three low- 
est modes (Cfc, w^) of each kind for the half-full tank. The order of magnitude 
of Ck and wk proves the weak intercoupling between shell and surface modes. 

The coupled frequencies of the full, three -quarter-full, half-full, 
quarter-full tank up to the seventh are compiled in Table I. To show once 
more the low coupling effect both kinds of frequencies are also determined 
without consideration of coupling and presented in Table I. The uncoupled 
liquid surface frequencies are obtained, assuming the tank is rigid, by the well 
known formula [ 8] . 
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FIGURE 2. LARGE ALUMINUM TANK, H/L = 0. 5; SECOND 
LIQUID SURFACE MODE, f 2 = 0. 59 
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FIGURE 3. LARGE ALUMINUM TANK, H/L = 0.5; THIRD 
LIQUID SURFACE MODE, f 3 = 0.70 




FIGURE 4. LARGE ALUMINUM TANK, H/L = 0. 5; FIRST 
SHELL DEFLECTION MODE, f = 7. 3 




FIGURE 5. LARGE ALUMINUM TANK, H/L = 0.5; SECOND 
SHELL DEFLECTION MODE, f = 16 













FIGURE 6. LARGE ALUMINUM TANK, H/L = 0. 5; THIRD 
SHELL DEFLECTION MODE, f_ = 21 


TABLE I. LARGE ALUMINUM TANK, COUPLED AND UNCOUPLED 
FREQUENCIES AT DIFFERENT LIQUID LEVELS 

R= 198.0 in. (502.92 cm) h = 0.22 in. (0.5588 cm) 

L = 488.5 in. ( 1240.79 cm) E = 1. 03 x 10 7 lb/in? (7101.85 x 10 7 N/m 2 ) 

p = 9. 359 x 10 -5 lb sec 2 / in. 4 (0. 001 x 10 4 kg/m 3 ) 

p = 2. 69 x 10 -4 lb sec 2 / in 4 ( 0. 288 x 10 4 kg/m 3 ) 
s 



H/L 

1. 


Number 
of mode 

Coupled 


1 

3.90 


2 

10. 1 

u 

>> 

3 

14. 3 

o 

§ 

4 

17.5 

& 

0 

u 

5 

20. 1 

XI 

bO 

6 

22.4 

*rH 

xi 

7 

24.5 


1 

0. 435 

0 

2 

0.589 

>> 

3 

0.709 

o 

« 

0 

4 

0.811 

& 

0 

£ 

5 

0.902 

O 

6 

0.984 

i — i 

7 

1.06 



0 . 

75 

Un- 


Un- 

:oupled 

Coupled 

coupled 

3.90 

5. 11 

5. 11 

10. 1 

12. 4 

12. 4 

14. 3 

17. 0 

17. 0 

17. 5 

20. 5 

J°. 5 

20. 1 

23. 5 

23.5 

22. 4 

26. 1 

26. 0 

24. 5 

28.4 

28.4 

0.435 

0.435 

0.435 

0.589 

0.589 

0.589 

0.709 

0.709 

0.709 

0.811 

0. 811 

0.811 

0.902 

0.902 

0.902 

0.984 

0.984 

0.984 

1.06 

1.06 

1 

1. 06 

l 


o 

Coupled 

7. 33 

16.0 
21. 3 

25.4 
28. 9 
32. 0 
34. 8 

0.435 

0.589 

0.709 

0.811 

0.902 

0.984 

1.06 


50 

Un- 

coupled 

7. 33 

16.0 

21. 3 

25.4 
28.9 

32.0 
34.7 

0. 435 
0. 589 
0.709 
0. 811 
0. 902 

0. 984 

1. 06 


0 

Coupled 

12.4 
23. 5 

30.5 

36.0 

40.6 

44.6 
48.2 

0. 431 
0. 589 
0.709 
0.811 
0. 902 

0. 984 

1. 06 


25 

Un- 

coupled 

12.4 

23.5 
30. 5 

36.0 

40.7 

44.7 
48.3 

0.431 
0. 589 
0.709 
0. 811 
0. 902 
0. 984 
1.06 
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co. = Jgit, tanh k. H 
1^1 1 

Ji (k. R) =0 


i = 1, 2, 3, . 


The uncoupled liquid-shell frequencies are determined using the simple 
analysis of Section V. For the coupling effect to become visible at all, the frequen- 
cies of Table I are given by three digitals. However, only in the low frequency 
range can one observe some slight effects. 

Finally, it should be mentioned that application of equation ( 54) shows 
that for the case under consideration the influence of the liquid compressibility 
is negligibly small. 


Brass Container 


Table n contains data as well as experimental and theoretical results. 
First and second frequencies at different filling heights are presented. Since 
application of equation ( 54) shows that the liquid compressibility can be neg- 
lected, the analysis of Section V has been applied. According to the condition 
(44) the first frequency for the liquid levels H>4 in. has been calculated by first 
order approximation using the first term of equation (43) . 


Steel Tube 


Data as well as experimental and theoretical results are compiled in 
Table m. The matter for discussion is the first frequency at different liquid 
levels. Since condition (44) is satisfied for all levels a first-order approxi- 
mation for calculation of the frequencies can be used (equation (43) ] . On the 
other hand, application of the condition (55) , 


h 

2R 


E 

K = °- 32 


shows that for the case in hand the liquid compressibility cannot be neglected. 
However, Table IH presents both kinds of frequencies, assuming the liquid 
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TABLE EL BRASS CONTAINER, EXPERIMENTAL 
AND THEORETICAL FREQUENCIES 


R = 4. 5 in. (11.43 cm) h = 0.005 in. (0.0127 cm) 

L= 11.2 in. (28.448 cm) E = 15. 3 x 10 e lb/in. 2 ( 105 483. 5 x 10 6 N/m 2 ) 

p = 9. 359 x 10~ 5 lb sec 2 / in. 4 (0.001 x 10 4 kg/m 3 ) 

p = 79. 55 x 10" 5 lb sec 2 / in. 4 (0. 851 x 10 4 kg/m 3 ) 
s 


H 

fl (Hz) 


f 2 (Hz) 


in. (cm) 

Test 

Analysis 

A(%) 

Test 

Analysis 

A(%) 

10.9 (27.686) 

217 

213.0 

1.8 

566 

548.3 

3. 1 

10.5 (26.670) 

226 

220.6 

2.4 

586 

563.4 

3.9 

9.5 (24. 13) 

250 

242.4 

3. 0 

627 

604.8 

3.5 

8.5 (21.59) 

276 

268.8 

2.6 

678 

652.5 

3.8 

7.5 (19. 05) 

309 

301.4 

2. 5 

736 

707.8 

3.8 

6.5 (16. 51) 

352 

342. 3 

2. 8 

811 

773.5 

4.5 

5.5 (13.97) 

404 

395.3 

2. 2 




4.5 (11.43) 

478 

465.7 

2.6 




3.5 (8.89) 

576 

563.4 

2. 2 




2.5 (6.35) 

739 

707.8 

4. 2 




1.5 (3.81) 

996 

955.5 

4. 1 

. 
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TABLE III. STEEL TUBE, EXPERIMENTAL 
AND THEORETICAL FREQUENCIES 


R = 1. 5 in. (3.81 cm) 

L = 16. 0 in. (40.64 cm) 
h = 0. 01 in. (0.0254 cm) 

E = 29 x 10 e lb/in. 2 (199 955 x 10 s N/m 2 ) 
p = 9. 359 x 10~ 5 lb sec 2 /in. 4 (0. 100 x 1 O' 1 kg/m 3 ) 
c = 0 . 56 3 x 10 5 in. /sec (1. 43 x 10 5 cm/sec) 



Experiment 

Analysis 



Liquid 

Assumed Compressible 

Liquid 

Assumed Incompressible 

H 

in. (cm) 

fi (Hz) 

f! (Hz) 

Deviation from 
Experiment 
(%) 

fl (Hz) 

Deviation from 
Experiment 
(%) 

16 (40.64) 

451 

436 

3. 3 

502 

11.3 

15 (38. 1) 

481 

465 

3. 3 

536 

11.4 

14 (35.56) 

518 

498 

3. 9 

574 

10.8 

13 (33. 02) 

553 

537 

2.9 

618 

11.8 

12 (30.48) 

594 

581 

2. 2 

670 

12.8 

11 (27.94) 

652 

634 

2.8 

730 

12.0 

10 (25.4) 

717 

698 

2.4 

803 

12.0 

9 (22.86) 

793 

775 

2. 3 

893 

12. 6 

8 (20. 32) 

894 

872 

2.5 

1004 

12.3 

7 (17.78) 

1003 

997 

0.6 

1148 

14.5 

6 (15.24) 

1176 

1163 

1.2 

1339 

13.9 

5 (12.7) 

1373 

1396 

1.7 

1607 

17. 0 

4 (10. 16) 

1655 

1744 

5. 4 

2009 

21.4 

3 (7.62) 

2133 

2326 

9.0 

2678 

25.6 
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compressible and also incompressible. In this way the significance of the com- 
pressibility for the case in question is demonstrated once more. 

SECTION III. EQUATIONS OF MOTION AND BOUNDARY CONDITIONS 


The cylindrical shell is referred to a cylindrical system of coordinates 
r, 6, z ( Fig. 7) . The shell displacements in directions r, 9, z shall be w, v, 
u. The surface elevation of the liquid is denoted by £. The following simplify- 
ing assumptions will be made: 

1. The analysis is restricted to the axisymmetric case ( v = 0) . 

2. The components of the shell inertial forces in direction u will be 
neglected. 

3. At the top, the shell is free and hence no axial membrane 
forces are acting. 

4. The liquid is considered inviscous and irrotational. 

Under these assumptions the equation governing the motion of the shell 
wall [9] is given by 
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FIGURE 7. PARTIALLY LIQUID-FILLED CYLINDRICAL CONTAINER; 
COORDINATE SYSTEM; DISPLACEMENTS 



H w(z,t) + ph 8 ^ t2 (Z,t) - =P(r,z,t) 0 — z ^ H 


(i) 


while the liquid velocity field can be represented by a gradient field 


v 


V cp (r ,z,t) 


0 ss r ^ R 
0 ^ z ^ H 


where the potential cp is a solution of 

V 2 (p (r, z,t) = 0 


or 


V 2 cp (r , z,t) 


i 9 2 _cp 

c 9t 2 


( 2 ) 

(3) 


depending on whether the liquid is assumed to be incompressible or compress- 
ible [2,3, 10J . 

Neglecting static and second-order terms the pressure equation is given 

by 

9 cp 

P (r,z,t) - Po = - p-jjj" (4) 

where p 0 is an external pressure independent from the potential cp, In the case 
of free vibrations p 0 can be assumed zero. 

The boundary conditions to be satisfied by cp are 

dcp (R,z,t) _ dw 

dr “ 8t 

9<P(r,0,t) = 

9z 


( 5 ) 

( 6 ) 


The free surface condition depends on whether or not the surface eleva- 
tion of the liquid is taken into account. In the first case the pressure on the 
undisturbed liquid surface equals--apart from a constant--the "weight" of the 
surface wave. In the second case the pressure is constant and can be assumed 
zero. Hence 
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pgt (r,t) 


(7) 


-P 


8 '(£_ 

at 



alternate condition 


( 8 ) 


If the surface elevation t, is taken into account an additional condition — 
equivalent to equation (5) — must be satisfied 


a# r,H,t) = 9£ 

dz at 


(9) 


SECTION IV. FREE VIBRATION ANALYSIS OF THE 
LIQUID-SHELL-FREE SURFACE SYSTEM 


The vibrations of this system are governed by equations (1) , (2) , (4) , 
through (7) and (9) . Using matrix notations, these equations can be written as 
follows : 


Eh 






' 




— — 

R 2 

0 


w(z,t) 


P h 
s 

0 

a 2 

w(z,t) 

a 

<P(. R,z,t) 

0 

Pg 


&(r,t) 


0 

0 

at 2 

&(r,t) 

- p at 

<P(r,H,t) 

_ 


vV( r, Z,t) 


a 

w(z,t) 


a^(R,z,t) 

9r 

at 

Ur,t) 


d<P(r, H,t) 

_ dz _ 


aV i dcp_ _aV 0 0 ^ r ^ R 

3r 2 r 3r 8z 2 0 ^ z ^ H 



a < j°(r,0,t) _ Q 
dz 


/ 


Assuming 


16 



( 11 ) 


<p( r,z,t) = H 2 o ) <P(r, z) cos cot 
w( z,t) = H w(z) sin cot 
C(r,t) = H Ur) sin cot 


( 12 ) 

(13) 


equations (13) reduce to 


where 


1 

0 


w(z) 

, ( 

i 0 


w(z) 


<P( R,z) 







+ -£ 


0 

Y_ 


J(r)_ 

W sh ( 

0 0 


C(r)_ 

6 



_ JL R 2 pg . . 

M R ’ Y Eh ’ 


P h 

s # 2 _ 

pR ; w sh R 2 p c 


E 


(14) 


(15) 


w(z) 


d<p(R,z) 

= H 

dr 

Ur) 

d<p ( r , H) 


dz 


d 2 <p i 8<p . 8 V „ 0 - r =5 R 

vV(r,z)=y+--+ 8z 2-0 o < z ^ H 


(16) 


(17) 


9y(r.O) =0 

8z 


(18) 


Equations ( 14) through ( 18) represent the eigenvalue problem in hand. 
The eigensolutions consist of the eigenvalues and the corresponding pairs of 
eigenfunctions : 


co 


k ’ 


w fc (z) 

4 k (r) 


k = 1,2,3. 


(19) 


which can be represented by the normal derivatives of a potential 


^ (k) (r,z) O^r^R; O^z^H 


( 20 ) 
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in accordance with equations ( 16) through (18) . On the other hand the potential 
<p(k) can be d er ived from the eigensolution ( 19) in accordance with the second 
boundary value problem of the potential theory. 


Let 


In the sequel, the orthogonality of the eigensolutions may be discussed. 


— - 


— «■ 

w 


w 

m 

; w 2 ; 

n 


n 






be two different eigensolutions. Then it follows from equation ( 14) in the usual 
way 

_6_ 


(co 2 - co 2 ) f w w dsj + co 2 if w dsj + f f dsA 

m n ^ m n m l n ^ n 2 I 


-co 2 (f w ^ (n) 
n V m 


dSi + / ? m ^ (n) dSz ) =0 


where 


2: r = R; 0 ^ 0 < 27T ; 0^z<H; dsj = Rd0 dz 

S: z = H; 0 < 0 < 2tt ; 0 — r — R ; ds 2 = rdr d0 


( 21 ) 


represents wetted shell wall and the free undisturbed liquid surface, respectively. 
Now , application of the well known Greens theorem shows that the two sums with- 
in the brackets are equal, hence the orthogonality relations can be written as: 


CO ^ CO 

m n 




( 22 ) 


Solving the eigenvalue problem will be performed using Galerkin's method. 
In doing so the potential <P will be approximated by 
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I 0 (kg r) cos kgZ 


with 


<P(r,z) 


n cos Kg z 

“ a f l^R) cos kgH 


n cosh Kg z >7 0 ( Kg r) 
+ ^ h i cosh K t H^o R) 


+ a 0 




( 23 ) 


J?7T 

H 


•^i^R) =0 


* = 1 , 2 , 


(24) 

(25) 


As it can be seen by substitution the series (23) satisfies Laplace’s 
equation (17) and the boundary condition (18). The three terms of the right 
side of equation (23) in succession represent the velocity potentials of: 


1. The liquid-filled container with flexible wall, rigid flat bottom, 
and rigid flat top on the undisturbed surface. 

2. The liquid-filled container with rigid wall and bottom and free 
liquid surface. 

3. The liquid-filled container with rigid bottom, free surface, and 

a ring -like flexibility of the wall. This potential governs the mean 
flow through the container wall and the undisturbed liquid surface. 


From equations (23) and (16) it follows 


<P( R,z) 

n 

cos kg z/ cos kg H 

n 

cosh ^ z/cosh Kg H 

<P( r,H) 

= Z 
£= 1 

I 0 (kgr)/I 0 (kgR) 

+ E b 4 
i= i 

)/y 0 (Kg R) 
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w(z) 


~R <P r (R,zf 

n 

rt 

cos k^z/cos k^H 


= M 


= M /j a 

Ur) 


R <F ( r , H) 
z 

1=1 

0 


( 27 ) 


n 

+ Z b , 


i=l 


1 


0 




+ a 0 p 


r i ^o(KgR) 


2p 


where 


cr( r) = r ; r( z) = z tanh z 

Io(r) 


(28) 


a £ = (T(k £ R) ;T I = t(k £ H) 5 i=1 * 2 > 3 - 


From equation (27) the pairs of coordinate functions of the Galerkin 
approach may be recognized as 


a. cos k.z/cos k.H 


0 


-1 

3 J J 

9 

T. 

• 

J 


0 


^ o(K.r)/J 0 ( K _R) 


2p 


j = 1,2,3,. . .n . 

To apply Galerkin 's method equations (26) and (27) must be substituted 
into equation (14). Now successive scalar (left) multiplication of equation (14) 
with the above pairs of functions and integration over wetted shell wall and liquid 
surface as defined by the domains (21) yields a linear homogeneous system of 
equations in the unknowns : 


au a 2 ; a 3 ;. . .a n ; bj; b 2 ; b 3 ;. . .b n ; a 0 . 


(29) 


In doing so, some of the integrations to be performed are the following: 



R 

J i o ( K.r) 3 o ( /c.r) ds 2 = 27 r J j7 0 ( ^ o ( K.r) rdr =. 

S 1 J 0 3 


(Orthogonality conditions) 
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cos k.z 

T L_ 

j cos k.H 
3 


cosh /Cg z 

cosh k. H 
l 



2(7. T. 

-yJ-i 
k. 2 + k 


RttH 


2 cos k.z cr. 

fa. r-L- dsj = 4jU 2 - J -g RttH 

J, j R cos k^H (k..H) 


r ^o(kt) Io(k.r) 

f — 1 —3 — Hso ~ 

J fx J 0 (k.R) M^R) 082 H 2 k| + k2 




RttH [11] 


Io(k £ r) °i 

/ Mk^R) dS2 = 4fl <k^ R ’ H 


mi 


t. / i r 2 \ ^o(K.r) 2 t 

/ ^ (■M 5i ) *»--? <^r *-> -tin 


In this way the problem will be reduced to the following symmetrical 
matrix eigenvalue problem 
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where 


Let 


a . = k.R 
] J 


j = 1,2,. . . 








<k) 


n 


a o 


(k) 


(31) 


be the eigensolution of the matrix eigenvalue problem (30). Then from equa- 
tions (23) and (27) the eigensolution (19) (20) follows as 


“k ; 



r n 

__ 

"( 

i=i 

n 

Z 

b i (: 

1=1 



<k) 003 V 

? ° l cos H a ° 


(k) 


) 


3 o ( ^ R) 


+ 2p 2 a 0 (k) 


_ $ (k) I » (k f r) 003 V $ (k> oosh 

‘ P k (T ' Z> ^ a < I 0 lk ( R) cos k,H + ^ 1 b i cosh R) 


+ a 0 


(k) 



0 < r ^ R ; 


0 < z < H . 
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The natural frequency may be obtained from the eigenvalue (31) as 



From equation (30) the orthogonality conditions of the eigenvectors (29) may 
be obtained easily. 


SECTION V. FREE VIBRATION ANALYSIS BASED ON THE 
ASSUMPTION OF ZERO ULLAGE PRESSURE 


The vibration is determined by equations (1) , (2) , (4) through (6) and 
(8) . Using the notation (11) , (12) and (15) these equations can be written as: 


w(z) =^T Jjv(z) 

(32) 

2 v B 2 (p 1 9<P 8 2 <p 

V cp(r,z) - gr 2 + r Sr + 8z 2 -0 

(33) 

W, Z ,=H 

(34) 

V(r,H) = 0 

(35) 

d<P(v , 0) 

9z U 

(36) 


Equations (32) through (36) represent the eigenvalue problem in hand. 
Solutions of equation (33) satisfying the boundary conditions (35) and (36) are 
given by 


<P( r,z) =I 0 (k £ 'r) cos k^’z 

£ = i,2,... 

_ (2 1 - 1)tt 
t 2H 


(37) 

(38) 
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From equations (32) (34) (37) the eigenvalues may be concluded as 

I 

Scr, 

,2 _ 


oo. = oo , — , — r 
l sh 1 + 6cr^ 


a l = ^(k £ R) 


(39) 

(40) 


where cr is defined by equations ( 28 ) . 

From equations (34) and (37) the correspondent mode follows as 

w.(z) = k^H Ij(k. ? R) cos kjz j = 1,2, — 

From equations (28) and (40) one recognizes 


lim t 

( J = oo 

i— o° l 

Hence equation (39) shows that the empty shell frequency ojg^ represents the 
upper limit of the frequencies oo ^ . 

In general the shell wall mass (as represented by p gh ) is small. Hence, 
for the lowest frequencies , 


6a ^ « 1 

and can therefore be neglected against 1. hi that case equation (39) reduces to 


"i = "sh 6 a l 


(41) 


Using the power series expansions of Io(r) and Ij(r) the first terms of 
the series expansion of o*(r) can easily be found. One obtains 


"< r > - T [‘ - i(§) ! + 1(f) 4 - « (l)‘ + • 


(42) 
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Although <x( r) is defined for all values of r the convergence of the series (42) is 
sure if r is at least smaller than or equal to 2. 


From equations (4i) and (42) it follows: 


k' 2 R 2 


00* = CO -T <5 
l sh 




(43) 


Formula ( 43) can be used for calculation of a >. if the shell wall inertia 
is negligible (low frequency range) and if k^'R is sufficiently smaller as two. 
If, in addition 


/k *R\ 2 

- - 

/ k i K \ i 

(2i -1)7 r R 

\ 2 / 16 

2 H_ 


a first-order approximation can be used. 


(44) 


SECTION VI. FREE VIBRATION ANALYSIS BASED ON THE 
ASSUMPTIONS OF ZERO ULLAGE PRESSURE 
AND LIQUID COMPRESSIBILITY 


Equations of motion and boundary conditions are given by equations ( i) , 
(3), (4) through (6) and (8). Using the notations (li), (12) and (15) these 
equations can be represented as 


w(z) = 


7 2./, / 


= -^T |w(z) +*r cp(R,z) J 

W sh L 6 J 

(45) 

’ dr 2 r 8r dz 2 c 2 

(46) 

- TT 8?>(R,z) 

w(z) " H — d^~ 

(47) 

7<r,H) = 0 

(48) 

d7(r,0) 

3z 0 

(49) 
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Solutions of equations (46) which satisfy the boundary conditions (48) 
and (49) are given by 


I 0 (rjkj 2 - co 2 /c 2 ) cos k» z 


i = 1,2, 


(50) 


where is given by equation (38) . This can be proved by substitution. 

The additional equations (45) and (47) can be satisfied by selecting 




CO 


sh 


6 <t(r V kf - co 2 /c 2 ) 

1 + 8 <r( R \/ k! 2 - co 2 /c 2 ) 


(51) 


where cr is defined by equation ( 28) . 

As in Section V it can be concluded that for the lowest frequencies , oo^ , 
it may be 


6 o (r / kj 2 - co 2 /c 2 ) « 1 . 
Hence, in that case equation (51) can be simplified for 


~ 1X1 sh 6 °'( R \/ k i " /° 2 ) ' (52) 

Now, if equation (44) is valid, that means low frequency range of cases with 
large H (for instance, liquid-filled pipes) then, in accordance with the series 
expansion (42) , one may replace equation (51) by the first-order approxi- 
mation 


co 


2 

i 
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Elimination of co. leads to 


a \ = k i 


i 2R \ 
K + Eh / 


(53) 


where K = pc 2 represents the bulk modulus of elasticity of the liquid and 



2R \ 
Eh / 


the pressure wave velocity of the liquid within the elastic container f 12] . 

Comparison of equations (39) and (57) shows that the compressibility 
of the liquid can be neglected if 




« 1 


C V 


(54) 


Application of this inequality to the cases covered by condition (44) means, in 
accordance with equation ( 52) , 


SECTION VII. FORCED VIBRATIONS, COMPARISON 
WITH BLEICH'S RESULT 

In the following a forced vibration analysis of partially liquid-filled 
cylindrical shell containers with flat rigid bottoms will be presented. The 
liquid is assumed to be compressible, the ullage pressure shall be zero. 

Let 


x = A sin t 


(56) 



be the excitation of the container. If the body-fixed coordinate system of Figure 7 
is used, the external pressure field created by the excitation (56) is given by 

p 0 = p (H - z) x (57) 

The equations of motion and boundary conditions are given by equations 
(1), (3) through (6) and (8) with p 0 given by equation (57). Using the nota- 
tions (11) and (12) these equations can be written as: 


where 


£2 2 £2 2 n r / z > 

w( Z ) - 2 W( Z )=- 2^ k-n(l- H 

sh sh L ' ' 

)J (58) 

2 d 2 (p 1 dcp d 2 <p S 2 2 _ 

^ ' ar 2 r 8r 8 Z 2 c 2 ^ 

(59) 

— , v a^(R,z) 

w(z) " H ' 9r — 

(60) 

?(r,H) = 0 

(61) 

8<P( r,0) 

dz 

(62) 

0<r^R, 0 < z ^ H 


II 

txl> 

(63) 


and 6, u, oj , are defined by equation (15). 
sh 

Solutions of equation (59) which satisfy the conditions (61) and (62) can 
be represented by 


where 


( 65 ) 


k* 2 = k! 2 
] J 


n 2 


and kj is given by equation (38) . Proof can be made by 

To satisfy the remaining equations (58) and (60) 
expansion in terms of the function set (64) 


substitution. 

the following series 


<p e ( r,z) = rj Yj c - Io( k * r ) cos k - z (66) 

j=1 J J J 


will be used. Substitution of this expression into equation (58) and subsequent 
integration over z from zero to H, taking into account the orthogonality 
conditions 


yields 


where : 



i,j = 1,2, . . . , 


» 2 Q 2 1 

1 6 ^h (k | H) 2 k*RI 1 (k*R) 






1 


, (67) 


cr*=cr(k^R) 1= 1,2,... 


( 68 ) 


and cr is defined by equations (28). 

Equations (65) , (67) and (68) show that resonance will occur if 

f = l,2,3,... , 


where is defined by equation (51) . 

From equation (37) , (64) and (65) one concludes that the liquid com- 
pressibility can be neglected if 
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J2 2 

c k 


7 « 1 . 


In that case equations (66) and (67) reduce to: 


<?’( r,z) =vYj c 2 I o( k / r ) cos k / z 


i=i 


C £ = 


_2_ fl 2 1 1 

6 o> s * (k^H^'RMk/R) ± 


u ) : 


( 69 ) 


where is given by equation (39) . 

The total pressure as indicated between the brackets of the right side of 
equation (58) follows from equation (66) as 


<P - V 




c 


i 


I°(k/r) cos 



(70) 


H. H. Bleich [2] derives the total pressure applying another method of 
solution. He uses a space-fixed system of coordinates with the origin at the 
liquid surface at rest. In the following, Bleich 's solution will be compared 
with that given above. To facilitate the comparison, Bleich's solution will 
be referred to a coordinate system having its origin at the center of the contain- 
er bottom at rest. Also sign rules and some notations of the paper in hand will 
be utilized. Bleich's potential is then given by 


7 ] jc 0 sin/3 0 (H - z) I 0 (/3 0 r) + 


OO 



(71) 
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Co = - 


2 Ij(/ 3 qR) 


IX 


3 


lo (fioR) - II (/ 3 oR)| 0o 00 R 2 COS 00 H 


(72) 


and 


2j t (0,R) 


c i = - 


Mp o (0 k B) +3l(0 £ R) ] 0 £ 0 i R 2 cosh 0 £ H 


(73) 


£ 0 is the only real root of 


S2 2 = 



gg(jSR) 

1 + 6cr(/3R) 


( 74 ) 


where pi, 6, a; , and <r are defined by equations (15) and (28). 
sn 

i0 t (0 £ > 0) i = 1,2,3,... 

represent the imaginary roots of equation (74) . Between the 0 and the 0 the 
following relations exist: 

0O 2 =0O ) 2 + -fr =0/ - ^2 1 = 1,2,... . (75) 

Equation ( 72) shows that the denominator of the right side vanishes if 


0o = 


(2,j - 1) 7T 

2H 



or 




Substitution of this expression into equation (74) and comparison with equation 
(51) shows that in this case 


and this represents resonance. 
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To prove that the expressions (70) and (71) represent one and the same 
quantity one first assumes the equality of both expressions then multiplies both 
series by 


r Io (/3 0 r) cos k!z 

^ j = 1 > 2 , 3 , . . . 
or 

i = 1,2,3,... 

rioOV) cos k j z 


respectively, and integration over the domains 

0=£zs=H; O^r^R, 

taking into account the orthogonality conditions of the trigonometric and Bessel 
functions, results in equations (72) and (73). 


George C. Marshall Space Flight Center 
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Huntsville, Alabama, March 21, 1967 
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